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Abstract 

Dark matter may be charged under dark electromagnetism with a dark photon that kinetically 
mixes with the Standard Model photon. In this framework, dark matter will collect at the center 
of the Earth and annihilate into dark photons, which may reach the surface of the Earth and 
decay into observable particles. We determine the resulting signal rates, including Sommerfeld 
enhancements, which play an important role in bringing the Earth’s dark matter population to 
their maximal, equilibrium value. For dark matter masses mx ~ 100 GeV — 10 TeV, dark photon 
masses mA' ~ MeV — GeV, and kinetic mixing parameters e ~ 10“® —10“”^, the resulting electrons, 
muons, photons, and hadrons that point back to the center of the Earth are a smoking-gun signal of 
dark matter that may be detected by a variety of experiments, including neutrino telescopes, such 
as IceGube, and space-based cosmic ray detectors, such as Fermi-LAT and AMS. We determine 
the signal rates and characteristics, and show that large and striking signals—such as parallel 
muon tracks—are possible in regions of the plane that are not probed by direct detection, 

accelerator experiments, or astrophysical observations. 
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I. INTRODUCTION 


Dark matter may live in a dark sector with its own forces. This possibility has some nice 
features. For example, the dark matter’s stability may be ensured not by some discrete parity 
imposed by hand, but simply by its being the lightest fermion in the dark sector. If the dark 
sector contains an Abelian gauge symmetry, dark electromagnetism, the dark photon and 
the Standard Model (SM) photon will generically mix kinetically. This mixing is of special 
interest because it is one of the few ways for a dark sector to interact with the known particles 
through a renormalizable interaction and it is non-decoupling: a particle charged under both 
dark and standard electromagnetism induces this interaction at loop-level, and the effect is 
not suppressed for very heavy particles. In this way, this is a prototype for simplihed dark 
matter models with light mediators. The idea of a separate sector with its own photon [1, 2] 
and the further possibility of kinetic mixing between these two photons [3, 4] were first 
explored long ago, and the myriad implications for dark matter detection have recently 
attracted widespread interest [5, 6]. 

In this framework, dark matter will collect in the center of the Earth and annihilate to 
dark photons XX ^ A'A. These dark photons may then travel to near the surface of the 
Earth and decay to SM particles, which may be detected in a variety of experiments, from 
under-ice/underwater/underground experiments, such as the current experiments IceCube, 
SuperK, and ANTARES, and future ones, such as KM3NeT, IceCube II, DUNE, and Hy- 
perK, to space-based cosmic ray detectors, such as the current experiments Fermi-LAT and 
AMS-02, and future ones, such as CALET, ISS-CREAM, and others. The resulting signals 
of electrons, muons, photons, and hadrons that point back to the center of the Earth are 
potentially striking signals of dark matter. 

The possibility of dark matter signals from the centers of large astrophysical bodies 
was first proposed and investigated many years ago [7-16], and there have been important 
advances for the particular case of the Earth in recent years [17-25]. Typically these signals 
rely on annihilation to neutrinos, resulting in single-particle signals with a continuum of 
energies. In contrast, dark photons decay into two charged particles, which may be seen 
at the same time in a single experiment, and the total energy of these charged particles is 
equal to the dark matter particle’s mass, producing potentially spectacular results. 

A schematic picture of this chain of events is given in Fig. 1. A number of processes must 
be evaluated to determine the resulting signal. For the specific case of dark photons, it is 
tempting to simplify the analysis by making a number of assumptions. For example, one may 
assume that the dark matter capture and annihilation processes have reached equilibrium in 
the Earth and that the capture cross section has some fixed value, such as the maximal value 
consistent with current direct detection bounds. Alternatively, the calculations simplify 
immensely for dark matter masses large compared to all relevant nuclear masses, mx rrixi 
or dark photon masses large compared to the characteristic momentum transfer so that 
the interaction is point-like. We show that none of these assumptions are valid in the 
regions of parameter space of greatest interest; the large mx approximation may lead to 
errors of an order of magnitude for mx ~ 100 GeV, and the large approximation may 
also lead to mis-estimates of factors of a few for very light ~ MeV. To accurately 
determine the sensitivity of experiments to probe the relevant parameter space, we carry 
out a general analysis, without making these simplifying assumptions. An early exploration 
of dark matter accumulation on the Earth mediated by massless dark photons is Ref. [4]. For 
previous work exploring the case of massive dark photons, see Ref. [26] for the case of dark 
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FIG. 1: Dark matter is captured by elastic XN —)• XN scattering off nuclei, collects in the 
center of the Earth, and annihilates to dark photons, XX —>■ A'A'. These dark photons then 
travel to near the surface of the Earth and decay to SM particles, which may be detected by a 
variety of experiments, including neutrino telescopes and space-based cosmic ray detectors. As 
an example, we show IceCube and various signatures there resulting from A' decays to electrons, 
muons, and hadrons. We discuss the possibility that double tracks (showers) may be resolved 
spatially (temporally) in the detector. 


matter capturing in the Earth and annihilating into neutrinos, Refs. [27-29] for early work 
on celestial body capture of dark matter annihilating into dark photons, and in particular 
Ref. [30] for a description of the general framework of annihilation to light mediators that 
highlights the specihc case of solar capture and gamma ray signatures, which were later 
searched for by the Fermi-LAT collaboration [31]. Finally, recent work has highlighted the 
effect of self-capture [32, 33] and boosted dark matter [34]. 

These results are timely for several reasons. Dark photons have attracted signihcant 
interest and are probed in many ways, including direct detection experiments, accelerator 
and beam dump experiments, and astrophysical observables [5, 6]. The signals we discuss 
are detectable for dark photon masses ~ MeV — GeV and mixing parameters e ~ 
10“® — 10“^, an interesting and large region of parameter space that includes territory that 
has not yet been probed. These values of niA' can also produce dark matter self-interactions 
that have been suggested to solve small-scale structure anomalies [35-39]. The range of £ 
values are naturally induced, for example, by degenerate bi-fundamentals in grand unihed 
theories [40]. It was recently pointed out that combining kinetic mixings of this size with 
the self-interacting models for small-scale structure can also explain the excess of gamma 
rays from the galactic center recently observed by Fermi-LAT [41]. 

At the same time, this work motivates a new class of searches for current indirect de¬ 
tection experiments to discover dark matter. At present there are a number of landmark 
experiments, including those mentioned above, that are transforming the held of indirect 
detection with high precision measurements and increasingly large statistics. In many cases, 
however, their sensitivities for dark matter searches are clouded by uncertainties in astro- 
physical backgrounds. The signals we highlight here come from a specihc direction (the 
center of the Earth), cannot be mimicked by astrophysics and, in many cases, are essentially 
background-free. As a result, the processes discussed here provide an opportunity for both 
current and future experiments to detect a smoking-gun signal of dark matter. 
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II. DARK PHOTONS 


We consider the simplest model of dark matter interacting throngh dark photons. The 
low-energy Lagrangian is 

L = 

+ "^f{i^-qfe4-mf)f + X{i^-gxf-mx)X , ( 1 ) 

/ 


where X is the Dirac fermion dark matter, and A and A' are the SM and dark sector 
gange helds with held strengths F and F' and hne-strncture constants a = e^/(dvr) and 
= fi'x/(4vr), respectively. The snm is over SM fermions / with SM electric charges g/. 
Dark electromagnetism is broken and the mass m^, is generated by some mechanism, snch 
as the Higgs or Stueckelberg mechanisms, which we assnme otherwise plays a negligible role 
in the signals discussed here. Note that the dark matter particles X are stabilized not by 
some ad hoc discrete parity symmetry or even by dark charge conservation (which is broken), 
but by Lorentz symmetry, since X is the lightest fermion in the dark sector. 

After diagonalizing the gauge kinetic and mass terms, the physical states are the usual 
massless photon A, which does not couple to X, and the dark photon A' with mass = 
— e^, which couples both to X and to SM fermions with charge £g/e, where £ = 
e/ Vl — e^. We take the independent parameters of the theory to be 

mx , ruA ', £ , ax • (2) 


We typically fix ax by requiring X to saturate the observed dark matter density through 
thermal freeze out, so ax = a^ ~ 0.035(mx/TeV). Alternatively, the maximum allowed 
coupling is set by bounds on distortions to the cosmic microwave background [42-44]. Fitting 
the results from Ref. [45], we find a™'^^ ~ 0.17(mx/TeV)^'®^. With a choice of ax the model 
is completely determined by the first 3 parameters. 

Dark photons decay to SM fermions with width 


r(A' ^ //) = 


NcS^qja{m\, + 2rnj) j^ 


3mA' 


m 


( 3 ) 


A' 


where Nc is the number of colors of fermion /. The dark photons we consider are produced 
from the annihilation of extremely non-relativistic X particles, and so have energy mx- For 
mA' S> me, the dark photon decay length is therefore 

-9X 2 


L — RmBe 


3.6x10 


/ mx/mA' \ f GeV 

V 1000 ) 


( 4 ) 


where ~ 6370 km is the radius of the Earth, and Be = B{A' —>■ is the branching 

fraction to electrons. The A' branching fractions can be determined from hadron production 
at colliders [46]. For mA' < 2m^, Be = 100%. As increases above 2m^, the 
A' —)■ /i+/i“ decay mode opens up rapidly, and Be drops to 50% at mA' ~ 300 MeV. For 
500 MeV < mA' < 3 GeV, Be and B^ are nearly identical and typically vary between 15% 
and 40%, with the rest made up by decays to hadrons, which also produce photons and 
neutrinos from meson decays. For mx at the weak-scale and mA' ~ 100 MeV — GeV, the 
requirement L ~ Rq, implies £ ~ 10“® — 10“^, and we will see that this is indeed the range 
of kinetic mixing parameters that gives the most promising signals. 
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III. DARK MATTER ACCUMULATION IN THE EARTH 


Dark matter interacting throngh dark photons is captnred and annihilates at the center 
of the Earth. The nnmber Nx of dark matter particles in the Earth obeys the equation 


dNx 

dt 


n _ r* 

^ cap 


ann-^^X 


( 5 ) 


where Ucap and Tann = ^C'ann-^x are the rates for the capture and annihilation processes. We 
ignore dark matter evaporation, which is negligible for weak-scale dark matter masses [11, 
12]. We also ignore self-capture from dark matter-dark matter self-interactions. The impact 
of self-capture for the Earth is suppressed by the fact that the escape velocity is low compared 
to typical galactic dark matter velocities, and so typical dark matter self-scatterings simply 
replace one captured dark matter particle with another [47]. 

The solution to Eq. (5) is 


T 


ann 


lce.ptatih" (^) 


( 6 ) 


where r© ~ 4.5 Gyr is the age of the Earth, and r = (GcapC'ann)”^'^^ is the timescale for the 
competing processes of capture and annihilation to reach equilibrium. To evaluate Tann, we 
must therefore evaluate both Gcap and Gann, which we now do in turn. 


A. Dark Matter Capture 


Dark matter particles are captured when elastic scattering off nuclei N in the Earth 
reduces their velocity below the escape velocity. The elastic scattering process XN ^ XN 
is mediated by f-channel A' exchange. The most relevant scattering targets, iV, are iron and 
nickel; these and other elements are listed in Table I. In the center-of-mass frame, the cross 
section is 


dcTjv 1 2e^ax(y.Z‘j^ 2 

dfl CM {Ex + [2p‘^{l - cos6cm)^ 

X [{ExEx +p^f + {ExEx + p'^ cosOcuf - {rux + m%)p'^{l - cos6'cm)] , (7) 


where Ex, Zx, rnx, and Ex are the energy, electric charge, mass, and nuclear form factor 
of target nucleus N, and p is the center-of-mass 3-momentum of the dark matter. Since 
the collision is non-relativistic, p is negligible everywhere in Eq. (7), except possibly the 
denominator. The cross section may then be simplified to 


d(T X 

dVt 


CM 


Ae'^axaZx 


(2p2(i 


^^N _ 

cos^cm) + 


Ex 


2 


( 8 ) 


where fix = + ^x) is the reduced mass of the X-N system. 

It is tempting to simplify the denominator by neglecting p, and reducing the A' exchange 
to a contact interaction. However, it is not always true that m\t ^ p^ so that the latter 
term may be neglected. The typical size of the momentum is p^ ~ fi'x'w'^, where w is the 
X velocity in the lab frame. Since capture typically occurs only for very small asymptotic 
dark matter velocities, a reasonable choice would be w = v^{rx) ~ 5 x 10“^, the escape 
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velocity at the radius that maximizes the radial number density nN{r)r‘^ of target nucleus 
N. With these values, the contact interaction limit fails for < 3 MeV. Rather than 
neglecting the momentum term altogether, a slightly more sophisticated approach would be 
to make the substitution p^(l — cos6*cm) —t In this work, however, we keep the full p 

dependence in the propagator and evaluate the capture rate numerically so that our results 
are valid throughout parameter space. We have conhrmed that our results reproduce those 
in the literature in the corners of parameter space where simplifying assumptions are valid. 
For example, they match Ref. [48] in the large-m^', point-like cross section limit. 

To determine capture rates, it is convenient to re-express the differential cross section 
as a function of recoil energy Er = — cos 6 'cm)/'^ 7 v in the lab frame. In the non- 

relativistic limit the expression simplihes to [49] 


d(7 ly 
dEji 




rriN 


w‘^{2m]yER m\,] 




N\ 


( 9 ) 


For Ejy, we adopt the Helm form factor [50], 


— exp [—Eji/Ex] , 


( 10 ) 


where Ejy = 0.114 GeV/H^^ is a characteristic energy scale for a target nucleus with atomic 
mass number A^- 

From this fundamental cross section we can determine the capture rate. The differential 
rate of dark matter particles scattering off nuclei with incident velocity w at radius r from 
the center of the Earth and imparting recoil energy between Er and Er + dE^ is given by 


da 


dGcap = nx V nN{r)^^^w /©(w, r) d'^w d-^r dEn , 


N 


dE 


R 


( 11 ) 


where rix = (0.3 GeV/cm^)/mx and nx{r) are the dark matter and target nucleus number 
densities, respectively, and f^{w,r) is the velocity distribution of incident dark matter at 
radius r, which is distorted from the free-space Maxwell-Boltzmann distribution, /(u), by 
the Earth’s motion and gravitational potential. We follow the velocity notation introduced 
by Gould [13, 14, 16] where n©(r) is the escape velocity at radius r and u is the dark matter 
velocity asymptotically far from the Earth. 

The total capture rate is obtained by integrating Eq. (11) over the region of parameter 
space where the hnal state dark matter particle has energy less than mxv^{r)/2 and is 
thus gravitationally captured. The escape velocity n©(r) and number densities nxir) are 
determined straightforwardly from the density data enumerated in the Preliminary Reference 
Earth Model [51]. Following Edsjo and Lundberg [20], the target number densities are 
modeled by dividing the Earth into two layers, the core and the mantle, with constant 
densities and elemental compositions given in Table I. The capture rate is then Gcap = 
^jyG^p, where the rate on target N is 


fRfi 


C'i = nx 


dr Allr'^nxir) / dwATTW^f^{w,r) 


dEn^ 0(AE) . (12) 


^dE 


R 


Here Q{AE) = Q{E m^^ — Eyr,\u) imposes the constraint that capture is kinematically possible 
by enforcing that the minimum energy transfer, Emin, to gravitationally capture the dark 
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Element 

Core ME 

Mantle ME 

Eilp(s-i) 

Element 

Core ME 

Mantle ME 

EilpCs-i) 

Iron 

0.855 

0.0626 

9.43 X 10^ 

Chromium 

0.009 

0.0026 

8.98 X 10® 

Nickel 

0.052 

0.00196 

7.10 X 10® 

Oxygen 

0 

0.440 

4.03 X 10® 

Silicon 

0.06 

0.210 

2.24 X 10® 

Sulfur 

0.019 

0.00025 

2.41 X 10® 

Magnesium 

0 

0.228 

1.05 X 10® 

Aluminum 

0 

0.0235 

1.62 X 10® 

Calcium 

0 

0.0253 

9.06 X 10® 

Phosphorus 

0.002 

0.00009 

2.04 X 10^ 


TABLE I: Mass fractions of the Earth’s core and mantle for the elements most relevant for dark 
matter capture [20, 52], Also shown for each element is the capture rate for mx = 1 TeV, 
mA' = 1 GeV, e = 10“®, and ax = — 0.035 as a measure of the relevance of the nuclear target 

for dark matter capture. 


matter particle is smaller than the maximum recoil energy kinematically allowed, i^max, 
corresponding to cos6*cm = — 1- Explicitly, these energies are 

1 2 ;/^ 

Emin = xmx [vo^ “ Emax = - -W^ . (13) 

2 Tflx 

To make further progress, we must determine the distribution fQ{w,r). A simple ap¬ 
proach is to only include the effect of the Earth’s gravitational potential. However, the 
Earth is within the gravitational influence of the Sun, and one might expect the acceleration 
of dark matter by the sun to suppress or eliminate the capture of heavy dark matter particles 
by the Earth. In 1991, however, Gould argued that the interactions of dark matter with 
other planets leads to diffusion of the dark matter population between bound and unbound 
orbits and one could thus ignore the impact of the Sun’s gravitational held and treat the 
Earth in the “free-space” approximation to reasonable accuracy [53]. 

More recently, however, this simple picture has been rehned with both potentially positive 
and negative implications. In numerical work, both Lundberg and Edsjo [20] and Peter [21- 
23] have investigated the inhuence of the Sun, Earth, Jupiter, and Venus in more detail, 
tracking the possibility that the Earth’s dark matter population is suppressed when particles 
are kicked out of the solar system or captured by the Sun. For the case of supersymmetric 
WIMPs—that is, dark matter with weak-scale mediators—they have found that these effects 
can reduce the Earth’s capture rate by one order of magnitude or more, depending on the 
dark matter mass. On the other hand, simulations of galaxies with baryons have shown that 
dark matter substructures may be pulled into the disk and create a signihcant and relatively 
cold enhancement of the local dark matter density known as a “dark disk” [54, 55]. For the 
case of WIMP dark matter, this population may enhance indirect detection signals from the 
Earth by up to three orders of magnitude [24, 56]. Note that the dark disk has a velocity 
relative to our solar system that is ~ 1/5 that of the ordinary dark matter halo [57]. It is 
thus plausible that the dark disk populates a region in phase space more amenable to Earth 
capture without signihcantly enhancing the direct detection rate. 

As we show below, the dark photon case differs signihcantly from WIMPs, because both 
the capture and annihilation rates are highly velocity dependent. One consequence of this is 
that r© is typically larger than r in Eq. (6), as opposed to the conventional wisdom that the 
Earth has not reached its WIMP capacity. It is therefore not possible to simply extrapolate 
the conclusions of WIMP studies to the present framework. In addition, as our analysis is 
valid for general dark matter and dark photon masses, inaccuracies in the particle physics 
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modeling are greatly reduced, and the astrophysical uncertainties from dark disk and other 
effects are very likely the dominant uncertainties entering the signal rate derivation. These 
astrophysical phenomena are therefore clearly interesting and important, but are beyond 
the scope of the present work. Here, we use the free-space approximation, not because it 
is the last word, but because it provides a simple “middle ground” estimate, with both 
suppressions and enhancements possible. 

With the free-space assumption, we proceed as follows. By energy conservation, w and u, 
the incident dark matter particle’s velocities in the Earth’s and galactic frame, respectively, 
are related by 

+ v'^{r) . (14) 

The capture rate for a general dtaj^/dE^ can then be rewritten as 


"Rs 


C'i = nx 


drAT:r‘^nx{r) / duA7iu^f(^(u 

Jo 


V? + u|(r) 


'0 JO “ 

Here /©(u) is dehned to be the angular-averaged and annual-averaged velocity distribution 
in the rest frame of the Earth [58], 


dER^Q{AE). (15) 


/©(«) = ^ 


dcos 6 dcos cf) f 


-1 


+ {Vq + V(s cos + 2 u{Vq + Videos (p) cos , (16) 


where V© ~ 220 km/s is the velocity of the Sun relative to the galactic center, and H© ~ 
29.8 km/s is the velocity of the Earth relative to the Sun. Many-body simulations and other 
considerations suggest a dark matter velocity distribution in the galactic rest frame of the 
form [59-64] 


f{u) = No 


exp 


-’gal 


U 


kuQ 


-i k 


- 1 


0(Wal - “) > 


(17) 


where A'^o is a normalization constant, Ugai is the escape velocity from the galaxy, and the 
parameters describing the distribution have typical values in the ranges [48, 65] 

220 km/s < Uq < 270 km/s 450 km/s < Ugai < 650 km/s 1.5 < A; < 3.5 . (18) 

We use the midpoint values of each of these, namely, uq = 245 km/s, Ugai = 550 km/s, and 
k = 2.5. The truncated Maxwell-Boltzmann distribution is recovered for A: = 0. 

Upon inserting Eq. (9), the dE^ integral in Eq. (15) evaluates to 


'■ dax 27Te‘^axaZjf 
dEuxTr- = -^ 




Xn 


■ -h Ei{-Xx) 


--min 


where we use the 


dEn w^mxEN 

in •- ‘^N 

substitution variable xx and exponential integral function [66], 


Xx = 


2mxER + m\i 

2mxEx 


Ei(j) = 




dt — . 
t 


(19) 


( 20 ) 


The total rate is Ccap = E Cc^p = ^2'K^e‘^axanxY.x ^n^'^nEx) ^exp c/^p, where 

rR« 


= 

^cap 


drr‘^nx{r) / duuf^{u)Q{Axx) -h Ei(— xat) 

Jo 


-.min 


( 21 ) 


—max 


The capture rates U^p for various nuclei iV at a representative point in parameter space 
shown in Table I. 


are 




















B. Dark Matter Annihilation 


Once a dark matter particle is captured by the Earth, it repeatedly re-scatters, drops 
to the center of the Earth, and eventually thermalizes with the surrounding matter. In 
the case of the Sun, the dark matter thermalizes within the age of the Sun for X-proton 
spin-independent scattering cross sections greater than 10“®^, 10“®°, and 10“^^ cm^ for 
mx = 100 GeV, 1 TeV, and 10 TeV, respectively [23]. Similar studies of Earth capture are 
not available. However, we will hnd that, for the range of parameters where an observable 
indirect signal is possible, the direct detection X-proton cross sections are at least ap ~ 
10“"^® cm^, corresponding to X-iron cross sections of crpe ~ ■Z'|g(mFe/'^p)^o'p ~ 10“"^^ cm^, 
many orders of magnitude larger than required for thermalization in the Sun. We therefore 
expect dark matter to be thermalized in the Earth to an excellent approximation. 

For thermalized dark matter, the annihilation rate parameter Gann is [48] 


Gann (^ann^) 


GxraxP® 

STca 


3/2 


( 22 ) 


where ~ 13 g/cm° and ~ 5700 K are the matter density and temperature at the center 
of the Earth, respectively, cxann is the cross section for XX —>■ H'H', and v is the relative 
velocity of the interacting particles, which is double the velocity of either interacting particle 
in the center-of-mass frame. 

The thermally-averaged cross section is 


— (c’"ann'i^)tree {^s) y 


(23) 


where 


(crann^)tree 


TTttx [1 — 

rrix [l — m\,/{2mx)]‘^ 


(24) 


is the tree-level cross section [67], and (Ss) is the thermal average of the S'-wave Sommerfeld 
enhancement factor. This Sommerfeld enhancement factor [68] has been determined with 
various degrees of rehnement. An analytic expression that includes the resonance behavior 
present for non-zero can be derived by approximating the Yukawa potential by the 
Hulthen potential [69-71]. The resulting Sommerfeld factor is 


Ss 


TT sinh(27rac) 

a cosh(27rac) — cos(27r\/c — a^G) ’ 


(25) 


where a = v/{2ax) and c = Qax'mx/ij’^rnA')■ The thermal average is, then. 



d?v 

(27rnQ)3/2 






(26) 


where z'o ^ 


C. Equilibrium Time Scales 

In Fig. 2 we present results for the equilibrium timescale r = (GcapGann)”^'^^ in 
the {rriA'yS) plane for mx = 100 GeV without Sommerfeld enhancement, and mx = 
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FIG. 2: Contours of constant rjT^, the equilibrium timescale in units of the Earth’s lifetime, in 
the {mA',£) plane for mx = 100 GeV without (top left) and with (top right) the Sommerfeld 
enhancement factor, mx = 1 TeV with the Sommerfeld factor (bottom left), and mx = 10 TeV 
with the Sommerfeld factor (bottom right). The dark sector fine-structure constant ax is set 
by requiring fix — 0.23. In the green shaded regions, the Earth’s lifetime is greater than the 
equilibrium timescale, r© > r. 


100 GeV, 1 TeV, 10 TeV with Sommerfeld enhancement. The dark coupling ax is hxed 
by the thermal relic density. For mx, the parametric dependence of r on £ and mA' 

enters dominantly through Gann and is r ~ ~ m\,/e. This can be seen in the baseline 

values of the contours in Fig. 2. The bumps in the contours reflect the resonance structure 
of the Sommerfeld enhancement factor (S' 5 ). 

In the shaded green (upper left) parts of the hgures, the Earth’s dark matter population 
has reached its maximal (equilibrium) value, and so the annihilation rate is essentially 
determined by the capture rate, with Fann ~ |Gcap- As one moves down and to the right, 
however, the equilibrium timescale grows, and the population is eventually not maximal. We 
will see that when the population is not at its maximal value, the signal quickly becomes 
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undetectable. 

The Sommerfeld enhancement plays an essential role in reducing the equilibrium timescale 
and making the signal detectable in large regions of the plane. For capture, the 

typical velocity that enters has an irreducible contribution from the gravitational potential 
that accelerates dark matter as it falls into the Earth. Capture interactions, therefore, occur 
at the typical escape velocity in the Earth’s core, Vesc ~ 5.0 x 10“®. However, after the dark 
matter particles are captured, they sink to the core, and come into thermal equilibrium 
with the normal matter. As a result, the population of dark matter particles at the center 
of the Earth is even colder, with relative velocities vq ~ 1.0 x 10“® [TeV/mx]^^^- In the 
mA' C(xmx limit, the Sommerfeld factor of Eq. (25) becomes 

_ 271 ax/v 

X _ Q-2-Kaxlv 



For thermal relics. So is therefore typically ~ 27iax/v ~ 10^ — 10®. Sommerfeld enhancement 
therefore reduces the equilibrium timescale by factors of ~ 100 for mx ~ 100 GeV, as can 
be seen in Fig. 2 by comparing the top right and top left panels. This reduction on r goes 
to ~ 1000 for mx ~ 10 TeV. The Sommerfeld factor therefore plays an essential role in 
boosting the current Earth’s dark matter population and the dark matter signal [26]. 


IV. SIGNAL RATES AND CHARACTERISTICS 


After dark matter accumulates in the center of the Earth and annihilates to dark photons, 
the dark photons propagate outwards with essentially no interactions with matter. The 
characteristic radius of the thermalized dark matter distribution in the Earth is [48] 


rx 


f 3T0 
\27iGNP®rnx 



150 km 



( 28 ) 


An observer at the surface of the Earth or in low Earth orbit therefore sees the majority of 
dark matter annihilations take place within 1.3°'\/TeV/mx of straight down. 

The dark photons are highly boosted with energy mx- In the decay A' —)■ //, the 
characteristic angle between the direction of a parent A' and its decay products in the 
Earth’s rest frame is 


9 ~ tan 


-1 


m\, — im'j 


m 


X 


1/2 


m\, — Am'j 


. , (29) 

assuming mA' rnx- Much larger opening angles are possible, however, as discussed in 
detail in the Appendix. 

The indirect detection signal is therefore two highly collimated leptons or jets that point 
back to the center of the Earth within a few degrees. As we will discuss, in some cases, the 
two leptons or jets may be simultaneously detected, and possibly even seen as two different 
particles, in contrast to the standard neutrino-based indirect detection signals, where there 
is only one primary particle. In any case, the signal of high-energy particles from the center 
of the Earth distinguishes the signal from all possible astrophysical backgrounds, potentially 
providing a smoking-gun signal of dark matter if the event rates are large enough. 

We now determine the event rates and characteristics for two classes of experiments: 
under-ice/underground/underwater detectors, represented by IceCube, and space-based ex¬ 
periments, represented by Fermi-LAT and AMS-02. 
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A. IceCube 


Dark photons may be detected in IceCnbe if they decay in IceCnbe or just below it. 
Decays A —)■ e^e~,qq will be seen as showers, and for m^/ > 300 MeV, typically 15% - 
40% of the decays will be to muons [46] and be seen as tracks. The number of dark photon 
decays that can be detected by IceCube is 

iVsig = 2 Tann ^decay T , (30) 

where the factor of 2 results from the fact that each annihilation produces two dark photons, 
Aes is the effective area of detector, 

edecay = - 6 -^^®+^)/^ (31) 

is the probability that the dark photon decays after traveling a distance between and 
+ D, where D is the effective depth of the detector, and T is the live time of the 
experiment. 

To very roughly estimate the detection rates for IceCube, we expect that for mx ~ 1 TeV, 
all dark photons that decay within the instrumented volume of IceCube are detected, and 
so we take Agg ^ 1 km^ and D ^ 1 km. For lighter dark matter, say, mx ~ 100 GeV, 
the decay products may be lost between the photomultiplier strings of IceCube. But these 
should be seen with high efficiency in DeepCore [72], the subset of IceCube with finer string 
spacings and lower threshold, and so we also present results for the instrumented volume of 
DeepCore, with A^g ^ 0.067 km^ and D ^ 0.55 km. 

In Fig. 3 we present the number of signal events for mx = 100 GeV, 1 TeV, and 10 TeV 
in the {mA',£) plane. The bumpy features and closed contours are real physical features 
resulting from Sommerfeld enhancement resonances. Also shown are the regions of parameter 
space disfavored by existing bounds on dark-photon-mediated XN ^ XN scattering from 
direct detection experiments, such as XENONIOO and LUX [73], and X-independent bounds 
on dark photons from beam dump experiments and supernovae [5, 74-77]. We use the 
recently updated supernova cooling bounds in Ref. [77] which differ from earlier calculations 
for emission of light particles by including multiple pion exchange processes. As a result, 
these cooling bounds are about an order of magnitude weaker than previous limits based on 
single pion exchange. 

We see that the indirect detection signal discussed here probes regions of parameter 
space that are so far inaccessible by other methods. As anticipated in Sec. II, the indirect 
detection signal is largest for £ ~ 10 “® — 10 “^, where the A' decay length is comparable 
to i? 0 . For mx ~ TeV and e ~ 10“®, for example, Wig ~ 10,000 events over 10 years 
are possible in regions of parameter space that are otherwise currently viable. IceCube has 
collected roughly 7 years of data already, and so detailed analyses will either exclude large 
new regions of the {mA',£) parameter space or discover dark matter. 

For mx ~ 100 GeV, the indirect and direct detection sensitivities are comparable for 
ax between a^x The indirect and direct detection sensitivities are shown in the 

conventional {mx, crxn) plane in Fig. 4, where axn is the spin-independent X-nucleon cross 
section: qi^A^axu = l^n^xT with T = Xe. The indirect detection signal is suppressed for 
both large axn (large e, dark photons decay too soon) and small axn (small e, dark matter 
capture is too slow and the captive population does not equilibrate). Of course, the large axn 
are already excluded by direct detection experiments. Focusing on the small axn region, for 
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FIG. 3: Red: IceCube event rates for T = 10 years live time in the {mA',^) plane for mx = 
100 GeV in DeepGore (top left) and IceGube (top right), mx = 1 TeV in IceCube (bottom left), 
and mx = 10 TeV in IceCube (bottom right). The dark sector hne-structure constant is set to 
the value which realizes fix — 0.23. Green: Single event reach for the maximal dark fine- 
structure constant allowed by cosmic microwave background distortion bounds [45]. Blue: 

current bounds from direct detection [73]. Gray: regions probed by beam dump and supernovae 
constraints [5, 74-77]. 

mx > 100 GeV, the indirect detection signals probe cross sections as much as three orders 
of magnitude below the current bounds from direct detection experiments, such as XENON 
and LUX. 

Since axn ~ a given axn corresponds to a larger value of £ when assuming the 

thermal a^x versus maximal Ox^^ dark sector coupling. For this reason, the dashed ax^^ 
curves on the left-hand plot in Fig. 4 are sometimes above the solid curves on the 
{mx,o'xn) plane. This is because when going from a^x fo the additional e reach 

gained in indirect detection experiments is less than that in direct detection experiments. 
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FIG. 4: Comparison of indirect and direct detection sensitivities in the {mx,crxn) plane for niA' = 
100 MeV (left) and 1 GeV (right). The direct detection bounds are from the LUX collaboration [78]. 
In this regime the interaction is effective point-like in contrast to the low niA' region [73, 79, 80] 
in Fig. 3, where the direct detection bounds become independent of for low niA'- Also shown 
is the ‘neutrino floor,’ where coherent neutrino scattering affects direct detection experiments [81]; 
the dashed line is an extrapolation. 


The reason for this is straightforward: the lower bound on the IceCube reach is set by the 
condition that the tanh^(r 0 /r) in Eq. ( 6 ) is ‘saturated’ near unity, i.e. that dark matter 
capture and annihilation are in equilibrium. This is why the lower contours of Fig. 3 display 
the same resonances as Fig. 2. Since this condition is set by the geometric mean of the 
capture and annihilation rates, it scales differently from direct detection experiments which 
has the same parametric dependence as the capture rate. 

For sufficiently small dark matter masses see, for example, the discussion 

below Eq. (2). This is seen in the right-hand plot of Fig. 4, where the Ns = 10 (dashed 
green) contour disappears in the low-mx bubble. 

The detector’s effective area and depth D are, of course, dependent on the energy 
and type of the dark photon decay products, and a more detailed study of detector response 
is required to estimate these more accurately. This is beyond the scope of the present work, 
but we note here some basic considerations. Muons with energies ~ 100 GeV — TeV lose 
energy primarily through ionization and travel a distance 


w = 



a 

a + (3Eth 


(32) 


before their energy drops below a threshold energy Eth, where p = 1 g/cm^, a ~ 
2.0 MeV cm^/g, and /3 ~ 4.2 x 10“® cm^/g [82]. For = 1 TeV and Eth = 50 GeV, 
on average muons travel a distance = 2.5 km. Dark photons that decay to muons a 
km or two below IceGube may therefore be detected in IceGube, and so the effective depth 
of IceGube is a bit larger than 1 km. For mx ~ 10 TeV, the effective depth is larger 
still, although less than a naive application of Eq. (32) would indicate, as such high energy 
muons lose energy primarily through radiative processes. At TeV energies, the experimental 
angular resolution for muon tracks is less than a degree, providing a powerful handle to re¬ 
duce background. For the case of showers from electrons or hadrons, the angular resolution 
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FIG. 5: Lab-frame muon velocity differences for A' —)• and {mx,rnA') = (100 GeV, 500 MeV) 

(left), (1 TeV, 1 GeV) (center), and (10 TeV, 1 GeV) (right). Distributions are normalized to 100 
events and different values of mA' are shown for comparison. The top axes measure the time delay 
between the two final states per km between the decay and observation positions. 


is worse, but still sufficient to identify showers as up-going to within tens of degrees. In 
addition, because dark photons decay completely to visible particles, contained events are 
mono-energetic, with the total energy equal to mx- The angle and energy distributions of 
tracks and showers are therefore completely different from astrophysical sources, and provide 
powerful handles for differentiating signal from background. 

The dark photon signal has two primaries, which could in principle be identihed as a 
smoking-gun signal for the dark sector. In Figs. 5 and 6 we show histograms of the velocity 
difference (time delay) and opening angle (track separation) of the two muons produced in 
a dark photon decay. Details of the distributions are presented in the Appendix. Parallel 
tracks have been considered previously in the context of slepton production from high energy 
neutrinos in Refs. [83, 84] and have recently been searched for by IceCube [85]. As a 
benchmark for IceCube reach, the parameters mx = 1 TeV, mA' = 1-3 GeV, and e = 2x 10“® 
gives an expected 40 muon events in 10 live years. The center panels of Figs. 5 and 6 then 
show that over ~ 2.5 km between the A' decay point and the maximal detection distance, 
one expects a few events with timing separation of ~ 0.03 ns and ~ 20 m track separation. 

The timing separation is below the IceCube Digital Optical Module timing resolution of 
~ 5 ns [86]. The track separations are less than the ~ 100 m separations probed by current 
analyses [85], but they are also much larger than the ~ 1 m separations from SM neutrino- 
induced charm production. These results motivate looking for parallel muon tracks with 
0(10 m) separations, which would be an unambiguous signal of physics beyond the SM, and 
a spectacular signal of dark photons and dark matter. One step in this direction is the pro¬ 
posed PINGU upgrade which would densely instrument a subset of the IceCube/DeepCore 
detector [87]. However, for the proposed dark photon search, this comes at a large cost 
in available volume and propagation distance. A possible alternative direction to improve 
sensitivity to these parallel muon signals is to increase the detector density of the IceTop 
surface array. 
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FIG. 6: Lab-frame muon opening angles (bottom axis)/track separation (top axis) for A' 
and the same mx values as Fig. 5. Different values of m^' are shown for comparison. 
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FIG. 7: 10 year signal event rates for a space-based detector in low Earth orbit, such as Fermi/AMS. 
Golors and bounds are the same as Fig. 3. 


B. Fermi-LAT/AMS-02 


The dark photon decay products may also be detected by space-based cosmic ray de¬ 
tectors, such as Fermi-LAT and AMS-02. Though these are far smaller than IceCube, the 
dark photon may decay anywhere between the Earth’s surface and the detector, providing 
a partially-compensating enhancement to the rate. For Fermi and AMS, we follow the for¬ 
malism described above, but now use = 1 m^. Both Fermi and AMS are in low Earth 
orbit, flying 550 and 400 km above the ground, respectively. We choose D = 550 km in 
Eq. (31). Note that, after the dark photon decays, the resulting charged particles are bent 
in the Earth’s magnetic held by an angle 


e 


„ TeV L B 

Q_5°__ 

p 550 km 0.5 G ’ 


(33) 


where p is the particle’s momentum, L is the distance it travels, and we have normalized the 
Earth’s magnetic held B to an average value at the surface of the Earth. For mx ^ TeV, this 
dehection is less than the dispersion from the dark matter’s spatial distribution at the center 
of the Earth given in Eq. (28), but for mx ~ 100 GeV, this dehection may be signihcant, 
and the signal may arrive at an angle as large as 5° relative to straight down. 

The resulting event rates for such space-based detectors are given in Fig. 7 for a live 
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time of 10 years. For mx ^ 100 GeV, the parameter space that can be probed largely 
overlaps with that already probed by direct detection, bnt Fermi and AMS may set bonnds 
complementary to the existing direct detection experiments. For higher masses one may 
probe nnexplored regions in the dark photon parameter space reaching m^/ ~ 0{1) GeV 
between e ~ 10“® and 10“^. As a benchmark, consider the parameters mx = 1 TeV, 
mA' = 830 MeV, and e = 3 x 10“^, for which one expects Nsig = 10 signal events in 10 
live years. The velocity difference and opening angle distribntions are shown in the center 
panels of Figs. 5 and 6. For a primary propagation distance of ~ 300 km, this yields timing 
separations of np to tens of nanoseconds and separations of np to a kilometer. We therefore 
do not expect to see both primary particles from dark photon decay in Fermi or AMS. Of 
conrse, this is still possible: the A' may decay near Fermi or AMS; secondary photons from 
hadronic hnal states my happen to have little transverse momentnm; or the A' may decay 
far from the detectors to two charged particles, which are both bent by the magnetic held 
into the detectors. Althongh possible, all of these are highly improbable, and two-particle 
events are a small fraction of the total nnmber of single-particle signal events. An alternative 
possibility is when there is a small splitting between and 2m/. In this case the decay 
prodncts have small transverse momentnm by Eq. (29), at the cost of a rednced branching 
ratio. 

Last, the nnmber of signal events N^g does not take into account experimental efficiencies 
associated with each apparatns. For example, we have assnmed that the volnme of the 
International Space Station between the Earth and AMS does not affect the dark photon 
primaries, and, fnrther, that the hadronic prodncts of the dark photons are detectable. A 
more complete analysis of the Fermi/AMS reach will reqnire more realistic modeling and 
different triggers. 


V. CONCLUSIONS 

We have presented a novel method to discover dark matter that interacts with the known 
particles throngh dark photons that kinetically mix with the SM photon. The dark matter 
is captnred by the Earth and thermalized in the Earth’s center, and then annihilates to 
dark photons. The dark photons then travel to near the snrface of the Earth and decay. 
We have determined the signal rates withont simplifying assnmptions concerning the dark 
matter and dark photon masses. In viable regions of the model parameter space, thonsands 
of snch dark photon decays are possible in IceCnbe, and smaller, bnt still detectable, signals 
in space-based detectors snch as Fermi and AMS are also possible. 

As with traditional indirect detection signals that rely on annihilation to nentrinos, the 
dark photon signal points back to the center of the Earth, differentiating it from astro- 
physical backgronnds. In contrast to the nentrino signal, however, the dark photon decays 
to two visible particles. The dark photon signal is therefore even more striking, as it is 
monoenergetic if fnlly contained. In addition, in principle both particles could be detected 
simultaneously yielding, for example, parallel muon tracks in IceGube with separations of 
~ 0(10 m). We have shown distributions of these separations for representative points in 
model parameter space. 

As discussed in Sec. Ill A, the leading uncertainty in the signal rate predictions is from the 
capture rate analysis. The escape velocity of the Earth is not large, and so this capture rate 
is subject to detailed modeling, including the effects of the Earth, Sun, Jupiter, and Venus. 
In addition, a cold “dark disk” population of dark matter may signihcantly enhance the 
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capture rates. The implications of these effects for WIMP dark matter have been considered 
in Refs. [20-24, 56]; it would be interesting to determine their effects on dark matter with 
dark photon-mediated interactions. 

In this study, we have assumed the dark matter X is a Dirac fermion and the mediator is 
a dark photon that mixes only with the SM photon, and so couples only to charged particles. 
It would be interesting to consider cases where X is a pseudo-Dirac fermion or a scalar, and 
cases where the dark photon mixes with the Z (and so couples to neutrinos, for example), 
or is replaced by a scalar (for which the dark matter may also be Majorana). Dark matter 
that collects and annihilates at the center of the Sun is also a promising source of decaying 
dark photons and will probe different regions of parameter space [88]. 

Finally, the experiments have been modeled very roughly here; detailed analyses, prefer¬ 
ably by the experimental collaborations themselves, are required to evaluate the accuracy of 
the signal rate estimates. However, our conclusion that there are viable regions of parameter 
space that predict thousands of signal events indicates that there are certainly regions of 
parameter space where the indirect detection signals discussed here are the most sensitive 
probes, surpassing direct detection detectors, beam dump experiments, and cosmological 
probes. The possibility of discovering signals of dark matter that, unlike so many other in¬ 
direct detection signals, are essentially free of difficult-to-quantify astrophysical backgrounds, 
provides a strong motivation for these searches. 
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Appendix: Decay Product Distributions 

We summarize analytic results for the kinematic distributions of the dark photon decay 
products, presenting the forward velocity difference between the two hnal states and the lab 
frame opening angle, which may be used to determine the time delay and track separation 
between these objects in a detector. For simplicity, we assume the dark photon decays 
isotropically in its rest frame. Angular correlations will modify our distributions, but will 
not change the ranges of time delay and track separation, which are our primary interest. 
With this approximation, in the center-of-mass frame, the dark photon decay products 
are evenly distributed in cos6*cm, where 6*cm is the angle between the dark photon boost 
direction and one of the decay products. The value of a kinematic quantity k for hxed model 
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FIG. 8: Velocity difference, Art, of the two particles produced in A' decay as a function of the 
center-of-mass frame angle cos 0 cm for representative values of a and b. 


parameters is a function k of cos 0 cm- The distribution of these values / is 

/w= E 


cos 6 »^^ 


k'(cos 0^m)I 


(A.l) 


Throughout this appendix we consider two-body decays A' —>■ ff and dehne 


a = 


2rr^ 

niA' 


h = 


rriA' 

nix 


(A.2) 


1. Velocity Distribution and Time Delay 


In the Earth’s rest frame, the forward velocities of the particles produced in A' decay are 

_ Vl - 6^ ± yi - COS0CM ,. o', 

1 ± Vl — cos 0CM ' 

where we use natural units c = 1. The difference of these velocities is 


Am = m+ — m_ = 


26^ a/I — cos 0CM 


26^ \/l — a? cos 0CM 


1 — (1 — IA){1 — a^) cos^ 0CM 1 — (1 — a^) cos^ 6, 


CM 


(A.4) 


where the last expression is valid for 6 ^ 1, the values we are most interested in. We plot 
Am(cos0cm) in Fig- 8. Observe that Aw scales like 6^ for small b; this is also seen in Fig. 5, 
where the mx = 1 TeV and mx = 10 TeV plots are related by a simple rescaling. Further, 
the distribution is fairly insensitive to a = 2m£/mA> for mA> ~ GeV and for i = e,fi. 

For a given Am, the (dimensionful) time delay between the two decay products for a 
decay that occurs a distance L from the detection point is 



CM- CM+ 


LAu L ^ 

-—Am , 

CU-Ua c 


(A-5) 


where we’ve taken the limit of large boost so that m± —)■ 1. 


19 






































FIG. 9: Earth-frame opening angle between the two particles produced in A' decay as a function 
of the center-of-mass frame angle cos 0 cm for representative values of a and b. 

2. Opening Angle and Track Separation 

In the Earth’s rest frame, the angles 9± of the decay products relative to the A' decay 
direction are 


tan0± = 


±6\/l — sin 0 cm 


\/l — lA ± a/1 — cos 0CM 
The opening angle between the two decay products is therefore 

26a/1 — a? sin 0 cm 


A0iab = tan 0+ — tan 0_ 


1 — (1 — a^) cos^ 0CM 


(A.6) 


(A.7) 


where the last expression is valid for 6 1. The scaling A0iab oc b can be seen in the center 

and right plots in Fig. 6. The maximal opening angle is 


A0 


max 

lab 


26a/1 — at cos 0 cm = 0 , 

b ^ /l-2a2 

- at cos 0CM = \ ^ : 

a V 1 — a 




(A.8) 


We plot A0iab(cos 0cm) in Fig. 9. For large a, the opening angle is maximized at cos 0cm = 
0, consistent with the intuition that the largest opening angle should correspond to fully 
transverse decays in the center-of-mass frame. But for small a, this intuition does not hold: 
the maximal opening angle occurs for cos 0 cm ~ 1 , where one particle is emitted “backwards” 
in the A' center-of-mass frame so that its forward velocity is signihcantly reduced, enlarging 
the opening angle. In most of the range of cos0cm; A0iab ~ 26, but the maximal opening 
angle A0J^^ ^b/a occurs for large cos 0 cm ~ 1 — \a^- 

Finally we show the correlation between A0iab and Am in Fig. 10. These plots identify 
where one may use the combination of the opening angle and time delay to discriminate the 
two hnal state particles. 
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FIG. 10: The correlation between the velocity difference and the lab-frame opening angle of the 
two particles produced in A' decay for representative values of a and b. 
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